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ABSTRACT 


Mean  horizontal  vlndapeed  profiles  within  and  above 
a  plantation  of  Japanese  larch  were  obtained.  A  log-profile 
analysis  of  above-vegetation  wlndspeeds  yielded  a  wide  range 
of  values  for  the  roughness  length  parameter  (zG)  and  the  zero 
plane  displacement  height  (D),  with  these  two  parameters  being 
highly  correlated  with  each  other.  The  computed  Eulerlan  space 
scale  of  turbulence  within  the  vegetation  showed  deeper 
penetration  of  large  eddies  after  needle  fall  and  during  high 
winds.  Power  spectra  showed  that  at  the  base  of  the  plantation 
most  of  the  variation  In  windspead  was  associated  with  gusts  of 
about  100  meters  wavelength.  Power  spectra  at  the  most  dense 
portion  of  the  plantation  canopy  showed  considerable  modification 
due  to  the  tree  spacing. 
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AIR  FLOW  AND  TURBULENCE  CHARACTERISTICS  IN  A 


JAPANESE  LARCH  PLANTATION 


L,  Hi  Allen,  Jf i 


INTRODUCTION 

Air  flow  within  and  Immediately  above  vegetation  couple*  the  land-mas* 
surface  source  or  sink  to  the  larger  weather  system  above.  Work  at  this 
location  in  the  past  has  been  concerned  with  the  source  and  sink 
distribution  of  water  vapor,  CO2,  sensible  heat,  and  momentum  within 
various  agricultural  crops,  using  both  energy  balance  and  momentum 
balance  approaches.  This  work  la  a  study  of  the  turbulence  characteristics 
within  a  tall  vegetated  surface. 

In  1964  a  study  was  begun  of  wind  flow  within  end  above  a  plantation  of 
Japanese  larch  near  Ithaca,  N,  Y,  The  sensor  system  consisted  of  cup 
anemometers  and  heated-thermocouple  anemometers  for  measuring  wind  flow, 
and  asnlrated  thermopiles  with  a  time  constant  of  about  one  minute  for 
measuring  temperature  profiles  above  the  vegetation. 

The  Japanese  larch  was  spaced  meters  by  3  maters  and  had  a  mean  height 
of  1040  cm.  There  was  quite  a  bit  of  individual  variation  in  heights  of 
the  tips  of  the  trees,  with  a  few  extending  200  cm  above  this  height. 

The  fetch  in  the  westerly  direction  was  about  150  meters. 

The  trees  were  sharply  tapered  at  the  tops,  with  the  vop  1/3  of  the  plant- 
air  layer  fairly  open,  the  middle  1/3  quite  dense,  and  the  lower  1/3 
consisting  mostly  of  dense  dead  branches.  Some  of  the  air  flow  patterns 
and  turbulence  characteristics  to  be  discussed  later  will  reflect  the 
effect  of  this  structure  of  the  vegetation. 

The  wind  data  were  taken  on  October  30  and  31,  after  partial  needle  fall, 
and  on  November  10  and  12,  after  total  needle  fall.  The  effects  of  the 
needles  on  the  flow  pattern  and  turbulence,  while  not  drastic,  do  show  up 
in  the  resulting  data. 

The  total  "yield"  of  needles  from  this  plantation  was  2200  pounds  per  acre, 
By  the  evening  of  October  31,  approximately  2/3  of  the  needles  had  fallen. 
However,  the  remaining  1/3  were  distributed  well  enough  so  that  the  drag 
characteristics  of  the  trees  should  not  have  been  changed  too  much  from 
that  of  a  fully  needled  tree.  Certainly  on  a  tree  such  as  larch,  every 
needle  cannot  function  Independently  as  a  drag  surface. 


The  computed  total  needle  surface  area  per  unit  ground  area  (LAI)  was 
9,36.  This  gives  an  average  leaf  area  density  of  about  0,013  cm2/cm3t 
assuming  equal  distribution  over  700  cm  of  depths  No  attempt  was  made 
to  evaluate  the  surface  area  of  the  branches  and  trunks 

Four  Beckman  &  Whitley  cup  anemometers  and  6  Thomthwaite  cup  anemometers 
were  used  above  the  vegetation  Six  Hastings-Raydist  heated-thermocouple 
anemometers  with  probe  types  N-7  or  N-7B  were  used  within  the  vegetation. 

The  cup  anemometer  data  were  read  from  banks  of  electromechanical  counters, 
with  sampling  durations  of  10  minutes  usually.  The  heated-thermocouple 
signals  were  recorded  on  an  analog  magnetic  tape  recorder.  Later  these 
data  were  converted  to  digital  form  and  stored  on  magnetic  tape  for  further 
processing. 

Table  1  gives  the  height  of  the  sensor  systems  used  in  this  study.  The 
wind  data  were  taken  at  heights  ranging  from  1569  cm  down  to  115  cm.  The 
juncture  of  the  cup  and  the  heated-thermocouple  anemometer  system  was  at 
1040  cm0 

WINDSPEED  PROFILES 

Thirteen  rims  were  made  with  both  cup  and  heated- thermocouple  anemometers. 
Several  other  runs  were  made  using  cup  anemometers  only.  Table  2  contains 
the  average  wlndspeed  data.  The  wlndspeed  data  obtained  at  a  height  of 
1040  cm  were  different  for  the  heated-thermocouple  anemometer  than  for  the 
cup  anemometer,  a  fact  commonly  noted  even  though  the  heated-thermocouple 
anemometer  had  been  carefully  calibrated.  Figure  1  shows  the  average 
normalized  wlndspeed  obtained  from  the  October  30  and  31  data.  The  wlndspeed 
profile  appears  nearly  linear  within  the  vegetation  in  the  lower  2/3  of  its 
depth.  Figure  2  shows  the  same  thing  for  November  12,  1964,  All  data  were 
normalized  with  respect  to  the  lowest  cup  anemometer.  These  figures  also 
show  a  tendency  for  the  wlndspeed  to  Increase  toward  the  bottom  of  the 
vegetation.  Tentatively,,  this  behavior  is  explained  on  the  basis  of  three 
factors,  all  of  which  are  dependent  upon  the  low  density  of  plant  parts  near 
the  ground.  First ,  the  fetch  for  this  site  was  inadequate,  but  it  was  the 
best  that  could  be  found  in  the  Ithaca  area.  This  feet  would  allow  gusts 
to  blow  through  the  bottom  of  the  plantation.  Secondly,  there  were  occasional 
holes  in  the  plantation  where  a  pressure  pulse  or  gust  of  wind  could  penetrate 
to  the  bottom  and  from  there,  with  most  of  the  effect  appearing  where  there 
would  be  less  obstruction  to  flow.  Thirdly,  the  slope  of  the  site  had  at 
this  point  at  least  a  3Z  downs lope  in  the  direction  of  the  prevailing  winds 
(west).  Anyhow,  these  profiles  indicate  that  downward  transfer  of  horizontal 
momentum  through  the  vegetation  cannot  account  for  all  of  the  flow  within  the 
vegetation,  an  assumption  that  has  frequently  been  used  with  short,  dense 
vegetation 

The  cup  anemometer  data  were  subjected  to  computer-programmed  log-profile 
law  analysis  according  to  a  procedure  by  Covey  (1963),  which  was  adapted 
from  a  personal  communication  from  Robinson  and  Tanner  No  buoyancy 
corrections  were  made  since  all  of  the  cup  anemometer  data  were  obtained 
near  the  top  of  the  vegetation,  and  because  of  the  shape  of  the  temperature 
profiles  obtained.  The  temperature  profiles  showed  a  maximum  value  about 
1  meter  above  the  top  of  the  vegetation,  with  decreasing  temperatures  above 
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and  below.  The  decrease  below  was  attributed  to  evaporation  and 
transpiration  from  plant  parts;  all  data  were  taken  after  rains  which 
left  trunks,  branches,  and  needles  visibly  wet.  However,  the  maxima 
value  of  the  temperature  profile  at  about  1  meter  above  the  mean  height 
of  the  vegetation  is  an  anomally  which  may  be  due  to  taller  trees 
upwind  or  to  net  upward  movement  of  air  due  to  slowing  down  of  winds 
within  the  vegetation. 

Typical  values  for  the  zero-plane  displacement  height  and  for  the 
roughness  length  parameter  were  635  cm  and  112  cm  respectively.  Taole 
3  shows  the  results  of  the  log-profile  analysis.  Figure  3  shows  the 
relationship  of  u*,  the  friction  velocity,  to  vindspeed  at  the  uppermost 
height  of  1569  cm.  The  slope  of  the  best  fit  relationship  passing 
through  the  zero  points  is  0.189  (when  the  same  units  are  used  on  both 
axes)  as  indicated  in  Figure  3. 

Earlier  work  by  Stoller  and  Lemon  (1963)  showed  5  one  changes  in  surface 
characteristics  of  flexible  vegetation,  such  as  alfalfa,  wheat,  and 
corn,  as  reflected  in  the  parameters  zero-plane  displacement  height  (D) 
and  roughness  length  parameter  (z  ),  with  increasing  windspeed.  One 
would  not  expect  very  much  changein  the  surface  drag  characteristics’ 
of  a  stlffer  vegetation  such  as  Japanese  larch  with  increasing  windspeed. 
Certainly  the  zero-plane  displacement  height  woulo  not  be  expected  to 
change  very  much  with  windspeed  due  to  a  "bending  over"  or  deforming  of 
the  plant  structure. 

Figure  4  shows  the  effect  of  changing  windspeed  on  z  .  The  fit  is  not 
too  good,  but  a  definite  trend  exists,  with  a  correlation  coefficient  of 
0.48.  Figure  5  shows  the  zero-plane  displacement  height  plotted  against 
windspeed.  The  correlation  here  is  better.  However,  part  of  the  goodness 
of  fit  of  these  data  may  be  due  to  the  coupling  of  the  higher  windspeed 
data  with  the  situation  where  all  of  the  needles  have  been  dropped  fra 
the  trees. 

Figure  6  shows  the  very  strong  correlation  of  z  with  displacement 
height,  with  r  —0.95.  Any  further  work  with  multiple  correlations  can 
only  Increase  r  by  a  negligible  amount.  Also,  Figure  7  shows  the  strong 
correlation  of  z  with  u*,  and  Figure  8  the  strong  correlation  of  u*  with 
D.  Analyses  of  variance  of  the  data  in  Figures  6,  7  and  8  all  indicate 
that  the  regression  line  is  significant  at  the  ,01  level. 

The  regression  lines  all  follow  the  pattern  found  by  Stoller 
and  Lemon  (1963)  for  corn. 

The  reason  for  z  and  displacement  height  varying  together  with  windspeed 
may  be  an  artifact  of  the  method  of  caputation  rather  than  having  any 
real  physical  significance.  In  the  log-profile  law. 
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if  *  defined  drag  coefficient  remains  anywhere  near  constant  with 
windspeed,  then  it  would  be  required  that 

u* 

—  -  constant, 

tm  9 

u 

which  in  turn  would  require  (a  +  D)/z0  to  remain  constant.  Hence,  any 
change  in  D  in  any  fitting  procedure  wouldjrequire  a  change  in  z0  to 
maintain  a  reasonably  constant  ratio,  u*/u.  However,  the  z0  vs,  u* 
and  u*  vs.  D  relationships  are  not  easily  explained  from  the  form  of 
the  log  profile  formula. 

The  log  profile  parameters  u*/u.^g,  z0,  and  D  were  compared  also  by 

dates  (October  30  and  31  vs.  November  12),  which  reflected  needle 
amount,  and  by  wind  direction.  The  first  two  parameters  showed  no 
significant  variation  due  to  either  date  or  wind  direction.  Displacement 
height  did  show  significance  at  the  ,05  F-level  for  both  date  and  wind 
direction.  However,  this  result  may  be  confounded  by  windspeed 
dependence  as  shown  in  Figure  5, 

The  various  relationships  indicated  in  Figures  3  through  8  were  olso 
investigated  by  regression  analyses  for  October  30  and  31,  November  10 
and  12,  November  10,  and  November  12  alone.  The  general  tendencies 
were  supported  in  most  cases.  Especially  did  the  u*,  zot  and  D 
relationships  maintain  a  high  level  of  significance. 

STATISTICAL  PROPERTIES  OF  TURBULENCE 

Figure  9  shows  how  the  computed  Euler lan  space  scale  of  turbulence, 

(using  the  heated-thermocouple  anemometers) 

rm 

L^  ■  u  i  R(t)  dt 

changed  within  the  vegetation.  The  values  are  quite  constant  within 
the  vegetation,  regardless  of  whether  windspeeds  were  high  or  low,  or 
whether  needles  were  present  or  missing.  There  was  a  slight  tendency 
for  values  to  be  higher  at  the  higher  windspeeds  and  with  no  needles 
present. 

The  biggest  change  occurred  at  the  725  cm  height.  The  tremendous  increase 
here  is  attributed  to  the  more  open  vegetation  after  needle-fall  and  to 
the  deeper  penetration  of  eddies  before  they  are  broken  up.  The  range  of 
the  eddy  scale  was  about  2  to  8  meters.  As  we  shall  see  from  power 
spectra  given  later  this  corresponds  to  a  shorter  eddy  scale  than  one  is 
usually  most  aware  of  under  gusty,  wl'idy  conditions, 

Turbulent  intensity  (coefficient  of  variation),  skevness,  kurtosls,  and 
cross-correlation  coefficients  with  respect  to  the  top-most  heated- 
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thermocouple  Anemometer  were  computed  for  each  run  and  then  averaged. 

Table  4  lists  the  average  values  at  the  six  anemometer  heights. 

In  all  cases,  the  distribution  of  windspeeds  are  skewed  toward  lower 
windspeed  values,  as  indicated  by  a  positive  skewness.  The  greatest 
skewness  shows  up  at  the  725-cm  anemometer,  which  is  in  the  most  dense 
part  of  the  vegetation. 

The  average  kurtosis  shows  a  tendency  toward  a  platykurtlc  distribution 
in  all  cases.  However,  the  windspeed  at  the  topmost  anemometer  shows 
only  a  small  deviation  from  that  expected  by  a  normal  distribution. 

Again,  the  windspeed  at  725  cm  shows  the  greatest  deviation.  It  has  a 
very  definite  platykurtlc  distribution.  At  this  level,  apparently 
the  anemometer  was  responding  well  to  large-scale  eddies,  but  then 
behaved  more  like  the  lower  anemometers  once  the  gust  passed  by. 

For  some  of  the  runs,  the  cross-correlation  coefficient  was  actually 
negative  for  the  bottom  two  anemometers.  This  would  indicate  that  the 
vertical  distance  is  large  enough  for  the  variations  in  windspeed  to 
tend  to  get  out  of  phase  with  respect  to  the  topmost  anemometer. 

Time-lagged  cross-correlation  coefficients  should  be  computed  to  find 
the  time  required  for  eddies  to  penetrate  to  given  levels. 

POWER  SPECTRA 

Power,  or  energy,  spectra  were  obtained  from  the  windspeed  data 
following  a  modified  procedure  of  Blackman  and  Tukey  (1958),  as  presented 
in  Pasquill's  ATMOSPHERIC  DIFFUSION  (1962),  In  order  to  shorten  computing 
time,  the  spectra  obtained  are  composite  spectra,  obtained  as  indicated 
by  Griffith,  Panofsky,  and  Van  der  Hoven  (1956)  and  by  Jones  (1957)  M.R.P. 
#1044.  The  point  at  which  the  two  sections  overlap  is  at  0,3  cps  in  the 
following  figures.  The  spectra  have  been  corrected  to  the  first  approxi¬ 
mation  for  linear  trends  in  the  data  based  on  the  treatment  by  Webb  (1955). 
The  correction  was  applied  to  the  autocovariance  function.  Likewise,  the 
spectra  were  corrected  for  sampling  duration  and  for  averaging  time  as 
outlined  by  Pasqulll  (1962),  Prewhitening  of  the  data  was  tried  and 
rejected.  It  tended,  in  these  cases,  to  vastly  overemphasize  the  magnitude 
of  the  spectrum  function  at  low  frequencies-  The  technique  presented  by 
Pasqulll  (1962)  was  the  one  tried-  Aliasing  of  the  data  was  no  problem 
because  the  filtering  action  of  the  thermocouple  anemometer  prevented 
frequencies  higher  than  that  allowed  by  the  sampling  time  from  entering. 

The  response  time  of  the  anemometers  was  1  second,  whereas  the  sampling 
time  was  0,16  second-  No  attempt  was  made  to  estimate  the  statistical 
reliability  of  the  various  points  in  these  spectra. 


Figure  10  shows  the  spectra  obtained  during  the  first  10-minute  run  on 
October  30,  1964,  The  values  plotted  are  nF(n)  normalized  on  the  basis 
of  the  maximum  value  of  nF(n),  as  a  function  of  frequency  in  cycles  per 
second  The  windspeed  at  the  uppermost  cup  anemometer  (at  1569  cm)  was 
371  cm/sec.  In  each  of  the  spectra  at  the  six  heights,  there  is 
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*  pronounced  low-frequency  peek  et  about  0.04  cps,  corresponding  to  a 
period  of  about  25  seconds.  This  peak  la  associated  with  gusts  or  eddies 
of  about  100  meters  wavelength. 

One  striking  feature  of  the  spectrum  at  a  height  of  115  cm  is  that  there 
is  very  little  contribution  to  the  variance  of  wlndspeeds  at  higher  fre¬ 
quencies,  This  means  that  at  the  floor  of  this  forest  there  is  less 
turbulence  on  a  small  scale  and  that  most  of  the  variation  in  horizontal 
air  flow  is  due  to  pressure  waves  associated  with  larger  scale  eddies. 

At  the  higher  frequencies  there  are  some  less  well-developed  peaks  in 
the  spectra  appearing  between  periods  of  3  to  7  seconds.  Since  the  mean 
wlndspeeds  at  these  levels  were  about  100  co/sec,  and  were  probably  lower 
during  periods  between  gusts,  these  peaks  are  probably  due  to  local  eddies 
created  by  individual  trees.  The  spacing  between  trees  would  be  about 
3  to  4  meters,  depending  on  whether  the  wind  direction  was  in  line  or  on 
the  diagonal  with  the  rows  of  trees, 

Two  other  examples  of  power  spectra  at  the  6  heights  are  in  the  next  two 
figures.  The  power  spectra  shown  in  Figure  11,  obtained  on  October  31,  1964 
at  1430  to  1440  EST,  shows  the  same  trends  mentioned  in  Figure  10,  The 
mean  windspeed  at  1569  cm  was  341  cm/sec.  The  period  between  large  gusts 
was  about  21  seconds,  giving  an  eddy  scale  length  of  more  than  70  meters. 

In  both  Figure  10  and  Figure  11  there  appears  to  be  a  kind  of  transition 
zone  at  a  height  of  725  cm,  where  the  density  of  plant  parts  is  greatest. 
More  of  the  fluctuations  are  due  to  higher  frequency  components,  having 
periods  of  the  order  of  3-seconds,  The  average  windspeed  was  .iHout  60  cm/ 
sec  in  Figure  11,  which  Indicates  an  eddy  scale  of  about  180  cm. 

The  third  figure  of  this  series  (Figure  12)  shows  spectra  under  conditions 
of  high  windspeed,  about  837  cm/sec  at  a  height  of  1569  cm.  The  low 
frequency  spectra  maxima  occurs  at  a  period  of  about  10  seconds  -  which 
yields  an  eddy  scale  length  of  above  80  meters  This  figure  also  shows 
very  little  contribution  at  very  low  frequencies  (less  than  0,02  cps). 

The  spectrum  at  a  height  of  115  cm  has  a  shape  that  would  be  associated 
with  an  exponentially  decreasing  autocorrelation  coefficient, 

CONCLUSIONS 

In  conclusion,  three  points  can  be  emphasized-  First,  the  application  of 
the  log-profile  law  to  these  data  leads  to  a  lot  of  variation  in  z0,  u*  and 
D,  most  of  which  appears  to  be  closely  correlated  each  with  the  other.  Even 
though  D  is  lower  at  higher  windspeed,  this  effect  may  be  related  also  to 
the  lack  of  needles  on  the  stand  when  wlndspeeds  were  high.  The  fact  that 
D  may  be  quite  a  bit  lower  was  also  indicated  by  Figure  9  which  showed  the 
tremendous  increase  of  the  Eulerian  space  scale  of  turbulence  at  the 
725-cm  height-  This  latter  fact  was  the  second  point  of  this  study. 

Thirdly,  the  power  spectra  types  show  a  pronounced  peak  corresponding  to 
wavelengths  approaching  100  meters.  Most  of  the  variation  of  windspeed  at 
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the  bottom  of  the  atand  (115  cm)  appears  to  be  associated  with  these 
large-scale  gusts,  for  reasons  mentioned  earlier  in  the  paper,,  In  the 
region  of  maximum  density  of  plant  parts,  the  spectra  show  Increased 
relative  contributions  from  small-scale  eddies  which  have  length  scales 
corresponding  to  the  tree  spacing  distances. 

The  type  of  turbulence  data  presented  here  are  currently  being  analyzed 
for  several  agricultural  crops  (oats,  corn,  soybeans  and  sunflower). 
Eventually  relationships  between  the  vegetation  height  and  structure  and 
the  turbulence  and  exchange  processes  will  be  obtained. 
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Anemometer  System 

Beckman  &  Hastings 

Whitley  Thorn thwalte  Heated-  Temperature 

cup  cup  thermocouple  system 


1569 

1340 

1040 

1510 

1492 

1260 

875 

1280 

1416 

1180 

725 

1130 

1340 

1100 

575 

1030 

1060 

345 

1040 

115 

CO 

CM 

NO 

CM 

NO 

00 

00 

00 

00 

00 

co 

NO 

CO 

cc 

•3 

*3 

ON 

CO 

<r 

ON 

co 

co 

O 

rx» 

CM 

co 

*3 

CO 

*3 

•3 

CO 

CO 

CO 

CM 

CO 

CM 

m 

o 

oc 

oo 

co 

NO 

NO 

o 

ON 

CO 

O 

<r 

ON 

ON 

NO 

00 

CM 

vO 

ON 

pH 

O' 

NO 

NO 

o 

<r 

o 

r>- 

NO 

CO 

co 

00 

CM 

ON 

fv 

o 

o 

<r 

ON 

co 

•3 

CO 

•3 

CO 

CO 

CO 

CM 

CO 

CM 

<r 

NO 

CO 

00 

NO 

m 

co 

co 

00 

CO 

H 

NO 

00 

co 

in 

NO 

ON 

CM 

<r 

co 

CO 

ON 

in 

ON 

•3 

CO 

CM 

ON 

NO 

<r 

CM 

r>. 

pH 

CO 

co 

00 

co 

CM 

NO 

co 

<*> 

CO 

CO 

•3 

CO 

CO 

to 

CM 

co 

CM 

•*3 

NO 

r*% 

r- 

NO 

m 

cm 

CO 

00 

o 

CM 

00 

*3* 

CM 

CO 

o 

•3 

ON 

c- 

O 

CO 

CM 

<r 

CM 

m 

o 

00 

H 

CM 

On 

O 

vO 

o 

CM 

m 

pH 

co 

m 

o 

o 

*3 

co 

to 

CO 

CO 

<3* 

CO 

CM 

CO 

CM 

co 

CM 

•3 

NO 

r- 

o 

m 

m 

cn 

NO 

00 

NO 

00 

NO 

o 

CO 

CM 

oo 

00 

-3 

CM 

ON 

CO 

00 

in 

CO 

o 

CO 

<r 

CO 

pH 

co 

00 

ON 

CM 

NO 

r>* 

o 

CM 

CO 

CM 

CO 

co 

co 

CM 

CM 

CM 

CM 

CM 

*3 

m 

NO 

r- 

o 

m 

m 

CM 

CM 

m 

CM 

CM 

ON 

O 

00 

pH 

00 

ON 

O 

co 

r>» 

■*«i 

CM 

pH 

NO 

r>» 

H 

NO 

CO 

NO 

00 

NO 

o 

co 

o 

o 

pH 

m 

pH 

CO 

r-v 

CM 

CM 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

*3" 

in 

O 

o 

NO 

in 

m 

<r 

O' 

i-4 

NO 

•3 

O 

O 

O 

ON 

CO 

O 

pH 

NO 

*3 

ON 

O 

o 

co 

r>» 

CM 

o 

CM 

NO 

CO 

•■3 

pH 

ro 

ON 

co 

pH 

ON 

CO 

NO 

co 

CO 

CM 

CM 

CM 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

pH 

co 

m 

in 

NO 

in 

O 

CO 

ON 

o 

CM 

r>. 

NO 

CM 

NO 

o 

o 

o 

UO 

in 

*3 

n* 

o 

ON 

m 

O 

CO 

o 

<r 

rJ 

CO 

ON 

CM 

00 

NO 

CO 

VO 

O 

co 

o 

o 

CM 

CM 

CM 

CM 

CO 

cm 

CM 

CM 

p-H 

CM 

pH 

CO 

>? 

m 

NO 

in 

*3 

<3 

*3 

CO 

ON 

NO 

ON 

NO 

00 

o 

ON 

o 

CM 

co 

•3 

pH 

lO 

pH 

O 

<r 

ON 

NO 

CO 

CO 

o 

CM 

ON 

o 

in 

in 

CO 

CM 

m 

o 

CM 

CM 

i-H 

CM 

CM 

CM 

CM 

CM 

pH 

CM 

pH 

co 

-3 

m 

m 

in 

<r 

*3 

O 

CO 

0 

o 

ON 

ON 

CM 

00 

*3 

o 

ON 

CO 

CM 

ro 

c 

o 

<r 

*3 

NO 

CO 

*3 

00 

J 

CM 

CO 

CO 

CO 

<r 

•pH 

a) 

«5 

o 

p-H 

p— H 

p-H 

CM 

p-H 

p-H 

CO 

CO 

CO 

n 

0) 

u  m 

a 

a 

0 

w 

u 

E  m 

T3 

c3 

C 

C 

Z) 

•pH 

w 

s: 

u 

4J  m 

o 

i  <r 

a 

.O 

TJ  co 

no  o  <r 

O'  N  W 


O  CO  NO 

no  r-*  m 


in  ^  co 
in  no  in 


n  n  o\ 
m  m  <r 


co  *3  in 
m  m  <r 


on  oo  o 

o  cm  o 


N  NO  00 
on  r**  m 


co  no  m 
no  in 


N  VO  N 

in  m  ^ 


mow 

in  no  m 


co 

CM 

co 

r>* 

•>3 

co 

cm 

CO 

NO 

CM 

NO 

**3 

co 

co 

CO 

CO 

CM 

CM 

<r 

CM 

pH 

<r 

NO 

CO 

CO 

O 

r- 

O' 

CM 

CO 

CO 

co 

pH 

pH 

CO 

•3 

r- 

ON 

pH 

r>» 

pH 

pH 

CN 

m 

•*3 

pH 

CM 

CM 

pH 

pH 

pH 

00 

NO 

CO 

CM 

00 

r- 

*3 

co 

ON 

CO 

CO 

CM 

pH 

pH 

pH 

pH 

pH 

pH 

CM 

*3 

vO 

•3 

pH 

CM 

pH 

H 

O 

CO 

CO 

pH 

CM 

CM 

CM 

pH 

H 

cm  co  <r  in  u  c 

r«.  r^.  r^. 


ro  O'  O  ^ 

n  r\  rs  c0  cc 


Table  3.  Log  profile  parameters  for  1964  Japanese  larch 


Date 

Run 

No. 

Mean  wind 
direction 
(degrees) 

u* 

(cm/sec) 

z 

0 

(cm) 

/D/ 

(cm) 

U1569 

(cm/sec) 

Oct.  30 

70 

292 

57.7 

57.1 

809 

371 

Oct.  30 

71 

292 

84.8 

128.8 

633 

418 

Oct.  30 

72 

292 

62.7 

81.3 

750 

362 

Oct.  31 

A 

345 

57.2 

39.2 

784 

436 

Oct.  31 

73 

334 

119.8 

206.7 

494 

487 

Oct.  31 

74 

334 

53.7 

52,2 

733 

372 

Oct.  31 

75 

299 

65.0 

112.4 

635 

341 

Nov.  10 

C 

307 

73.5 

169.0 

460 

341 

Nov.  10 

D 

292 

40.4 

42.9 

760 

296 

Nov.  10 

E 

292 

35.3 

15.1 

877 

338 

Nov.  10 

76 

319 

20.8 

5,8 

870 

248 

Nov.  12 

F 

193 

80.9 

74,1 

618 

517 

Nov.  12 

77 

200 

118,8 

96.6 

562 

698 

Nov.  12 

78 

200 

197,7 

221.6 

362 

837 

Nov.  12 

79 

200 

109.9 

42.3 

682 

838 

Nov.  12 

80 

200 

192.4 

223.3 

378 

807 

Nov.  12 

81 

200 

111.3 

83.7 

618 

678 

Nov.  12 

82 

200 

164,2 

274.7 

310 

628 

Table  4.  Averages  of  turbulent  intensity,  skewness,  kurtosis,  and 
cross-correlation  coefficient  with  respect  to  the  topmost 
heated-thermocouple  anemometer  at  the  indicated  heights  in 
Japanese  larch  vegetation,  1964, 


Height 

(cm) 

Turbulent 

intensity 

Skewness 

Kurtosis 

Cross¬ 

correlation 

coefficient 

1040 

=  47 

,50 

3,05 

1.00 

875 

,54 

1.11 

4,71 

.62 

725 

,57 

1,62 

8.16 

,44 

575 

c  51 

,94 

4,42 

.22 

345 

,57 

82 

4,47 

,19 

115 

-51 

.81 

4,92 

,13 
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AVERAGE  NORMALIZED 
WINDSPEED  IN  LARCH 
1964,  OCT.  30  -  31 


1,  Average  normalized  windspeed  above  and  within 
Japanese  larch  near  Ithaca,  J.  Y.,  Oct.  30  and 
31,  1964. 


Fig,  3.  Friction  velocity,  u*,  as  a  function  of  windspeed 
at  the  1569-cm  height  in  1964  in  Japanese  larch  near 
Ithaca,  N,  Y,  The  solid  dots  indicate  data  on  Oct. 30 
and  31,  the  crossed  open  circles  data  on  Nov.  10,  and 
the  open  circles  data  on  Nov,  12. 


0 
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JAPANESE  LARCH,  1964 


4,  Roughness  length  parameter,  z0,  as  a  function 
of  windspeed  at  the  1569-cm  height  in  1964  in 
Japanese  larch  near  Ithaca,  N.  Y,  The  solid  dots 
indicate  data  on  Oct.  30  and  31,  the  crossed  open 
circles  data  on  Nov,  10,  and  the  open  circles  data 
on  Nov.  12, 


DISPLACEMENT  HT.  (M) 


WINDSPEED1569  (M/SEC) 


Fig.  5.  Zero-plane  displacement  height  as  a  function 
of  wlndspeed  at  the  1569-cm  height  In  1964  in 
Japanese  larch  near  Ithaca,  N,  Y.  The  solid  dots 
indicate  data  on  Oct,  30  and  31,  the  crossed  open 
circles  data  on  Nov,  10,  and  the  open  circles  data 
on  Nov,  12. 
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DISPLACEMENT  HT.  (M) 


6.  Relationship  between  the  roughness  length  paramete 
zot  and  zero-plane  displacement  height  in  1964  in 
Japanese  larch  near  Ithaca,  N.  Y.  The  solid  dots 
indicate  data  on  Oct,  30  and  31,  the  crossed  open 
circles  data  on  Nov,  10,  and  the  open  circles  data  on 
Nov.  12. 


(CM) 


JAPANESE  LARCH,  1964 


Fig.  7,  Relationship  between  the  roughness  length  parameter, 
zQ,  and  the  friction  velocity,  u*,  in  1964  Japanese 
larch  near  Ithaca,  N.  Y.  The  solid  dots  indicate  data 
on  Oct.  30  and  31,  the  crossed  open  circles  data  on 
Nov,  10,  and  the  open  circles  data  on  Nov.  12. 
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0 


Fig,  8.  Relationship  between  the  friction  velocity,  u*, 
and  the  zero-plane  displacement  height  in  1964 
Japanese  larch  near  Ithaca,  N,  Y.  The  solid  dots 
indicate  data  on  Oct,  30  and  31,  the  crossed  open 
circles  data  on  Nov,  10,  and  the  open  circles  data 
on  Nov,  12. 
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POWER  SPECTRA  (nF(n)/ nF(n)max)  AT  VARIOUS  HEIGHTS  IN  JAPANESE  LARCH 
RUN  NO.  64-70  OCT.  30,  1964  1528-1538  EST 


.03  .06  .1  .3  .6  1.0  .03  .06  .1  .3  .6  1.0 

n  (  c  p  s  ) 


Fig.  10.  Composite,  normalized  power  spectra  at  the  indicated 
heights  obtained  from  heated-thermocouple  anemometer  data 
in  1964  in  Japanese  larch  near  Ithaca,  N,  Y. 


E 


n(c.p.s.) 


11,  Composite,  normalized  power  spectra  at  the 
indicated  heights  obtained  from  heated-thermo- 
couple  anemometer  data  in  1964  in  Japanese  larch 
near  Ithaca,  N,  Y, 


.03  .06  .10 
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03  oe  P 


q3  cpp 

n(cos; 


12.  Composite,  normalized  power  spectra  at  the 
indicated  heights  obtained  from  heated-thermo¬ 
couple  anemometer  data  in  1964  in  Japanese 
larch  near  Ithaca,  N.  Y, 


APPENDIX 


This  appendix  gives  a  listing  of  Lue  FORTRAN  computer  program  used 
to  compute  power  spectra  with  the  Cornell  University  Control  Data 
Corporation  1604  digital  computers 

This  program  was  designed  with  several  features  which  should  be 
pointed  outc  First,  it  was  set  up  with  the  ability  to  compute 
spectra  across  up  to  three  wavelength  ranges,  or  stages,  as  indicated 
in  the  program?  This  feature  reduces  computing  time  tremendously, 
and  prevents  overtaxing  of  computer  storage  capability* 

The  second  feature  permits  a  correction  to  be  applied  to  remove  the 
effect  of  a  linear  trend  in  the  data  upon  the  power  spectra*  The 
third  feature  allows  one  to  vary  a  "prewhitening"  factor  by  changing 
one  card*  The  fourth  feature  is  the  correction  for  finite  averaging 
time  and  finite  sampling  duration. 

This  program  was  designed  specifically  for  working  with  wlndspeed 
stored  on  data  magnetic  tape  in  a  particular  "format,"  However,  since 
it  is  written  with  many  subroutines,  it  is  possible  to  modify  the 
program  for  other  forms  of  data  input  and  other  types  of  data. 

The  subroutine  "CRACK"  was  written  by  David  Bessel,  Department  of 
Computer  Science,  Cornell  University,  Ithaca,  N,  Y, 

The  program  was  written  by  Al\in  Kaskel,  Therm  Advanced  Research, 

Inc. ,  Ithaca,  N,  Y 
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PROGRAM  PSPECTRA 


APPENDIX 


CONTROL  I OENT 1 F t  CAT  ION 


NSTA6E  ■  INTEGER  REPRESENTING  THE  PFRIOO  OF  THE 
NON-OVERLAPPING  means  TO  BE  analyzed 
NSTAGE* I  IS  TMF  RAW  DATA 

NSTAGECO  ■  INTEGER  EQUAL  TO  AN  NSTAGE  TO  WHICH  THE 

ANALYSIS  IS  TO  BE  COMPARED.  USED  to  DETERMINE 
THE  NUMBER  OF  POINTS  IN  THE  ANALYSIS 

INOExO  ■  I  USE  NON— ADJUSTED  COVARIANCES 
■  2  USE  ADJUSTED  COVARIANCES 

1CROPCHG  »  I  OATA  IS  FOP  SAME  CROP  AS  PREVIOUS  ANALYSIS 
■  2  OATA  IS  FOP  NEW  CROP 


PROGRAM  ALWAYS  ASSUMES  THAT  THE  STARTING  POINT  NUMBER  FOR  A 
RUN  IS  FOR  OATA  IN  CHANNEL  1  AND  THAT  THE  FINISHING  POINT 
NUMBER  IS  FOR  OATA  |N  CHANNEL  7 

type  Real  lh.nfn.nffn.nfnr.nffnr.npfn.npffn 

COKMON/NUMBERS/NBS.NHS.NBF.NPF .M.NP, NSTAGE. lNOexQ.NU.MOATA.TFAC. 

*  ofltat 

COMMON/CHANNELS/ 1 CHUSE  «  7 ) • l ANUSF 171. I  SAVE (71 

COMMON/XSUKS/XXOI.7)  «X| (31 .7).XK(3I  .  7  >  .  XS  (  1  000 . 7  >  •  DATA  (62 . 7  >  . 

*  XT (31. 7) .XJ(3I.7).TJ(3I .7) . XTBAR (7  > . S ( 7 ) 

COMMON/ ANEMFACS/U (233. 7) .ALPHA (7) .BETA (7) .A(7).B(7).C(7).0(7>. 

*  E(7).F(7) «G(7j«HC7) 

COMMON/VALUES/OK (31 ),0KK(31 ).LH(3I ).UH(3I ).FN0(3I ).FN(3I ).FFN(3l >• 

*  NFN(3I  ).NFFN(3l  ».NFNR(3I  I.NFFNRI3I  I.NPFNOl  ), 

*  NPFFN(31 ) . VH (31 ).FRE0(3I I.SIG(3l ).SuML 
COMMON/EXTRAS/VV (7) .VEL (7 | .VELS(7) .NSTAGECO 

COMMON/ I /NUMPER.NUMBLOOC.N'JMPNTS. IFLAG. ICHANNEL (1680 ) • IDATA (  1680) 
DIMENSION  1CHT(7) . 1WHT(7 J .FZ(7). IAN0(7> . ITCN0(7). I TNO ( 7 ) . 

*  CROP I OEN( 10) 

10CO  FORMAT (1613* 

|OlO  FORMAT (2F10.3.F 13.3. 3E 13.5/20X.AE13.3 > 

1020  FORMAT (315.F10.5) 

1030  FORMAT(IOAB) 

|OaO  FORMAT (BF | 0.3 1 

1030  FORMAT ( 1H1 .AOX.39HPOWER  SPECTRA  ANALYSIS  OF  RAW  WlNO  DATA) 

1 060  FORMAT( IHI .26X.46HP0WER  SPECTRA  ANALYSIS  OF  RAW  WIND  OATA  USING  • 

*  20H AO JUSTED  COVARIANCES! 

|070  FORMAT! 1H| . 34X . A AHPfJwER  SPECTRA  ANALYSIS  OF  PRE-WHITENEO  WIND  . 

*  7H0ATA.B* «F4.2 ) 

1080  FORMAT ( 1H1 .20X . 44HP0WER  SPECTRA  ANALYSIS  OF  PRE-WHITENEO  WlNO  . 

*  34H0ATA  USING  ADJUSTED  COVAP I AnCES «B« . F4 .2 > 

I  090  FORMATdH  .  1 9x .  1  0A8  .  /  ) 

|  1 00  FORMATdH  .2X.8HCR0P  HT. ,3X.1HZ.3X.4HF(Z> .3X.9HANEM.  N0..2X. 

*  8HTAPE  NO. .2X. I4HTAPE  CHAN.  NO . . 3X . 4HL AGS . 3X . 

*  I3HSAMPL1NG  T I ME . 3 X . 8HRUN  T |ME . 3X . 1 6HAVERAG I NG  PERIOD) 

I  I  I  0  FORMAT (  1 H  .I7.4H  CM..I6.4H  CM. .F6.4 • I  9. I  1 0 . 1 1  2 . 1 1  3 .FI 0.2 .5H  SEC.. 

*  F | C , 2 , 5H  SEC..F9.2.3H  SEC.) 

1120  FORMATdHO.4 1X.8HADJUSTED./. 14X.5HSIGMA.8X.7HQ  SUB  K.9X.7HQ  SUB  K, 

*  9X.7HL  SUB  H.9X.7HU  SUB  H> 

I  1 30  FORMAT! 1H  .3X.3F16.7) 
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1140  FOAM AT ( ) HO  >4 1 X  «  8H ADJUSTED. /< I4X  »9MS 1 GMA  ,8X «  7H0  SUB  F.9X.7M0  SUB  K. 

4  9X.7HL  SUB  H.9X.7HV  SUB  H.9X.7HU  SUB  H) 

1190  FORMAT! !H  .3X.6F16.7l 

1160  FORMAT! 1H0, 13X.9HSIGMA.6X. I IHFRE0UENCY.N.4X. 12HFIN)*0  SUB  0.10X. 

4  4MF IN  )  . 1 IX.6HN4FINI.4X. 16HN4F  «N>/IN*FIN)  )MAX >2X  « 

4  1 IHN449/34FIN) I 

1 170  FORMAT! 1H  .3X.7FI6.71 

1  ISO  FORMAT  1 1 HO • 99X • BH AO JUST ED<BX<6HADJUSTED<BX<8HADJUSTED<8X  . 

4  6H  ADJUSTED./  .  1 4X  •  9HS  1  CMA  •  6X  .11 HFREQUENC Y  <  N. 26X  .  4HF  INl.llX. 

4  6HN4FIN1.4X.IBHN4FIN1/IN4FIN)  IMAX.2X.  1  |HN**9/3*F  IN)  1 

1190  FORMAT! 1H  . 3X .2F 16. 7. 16X .4F 16. 7 1 

1200  FORMAT! IH0.97X.9HC0RR£CtED*/«9X< 1 IHSUM  L  SUB  M.6X.10HU  ICM/SEC). 

•  4X.14HSUM  OF  SQUARES. 2X.14HSUM  OF  SQUARES. 3X . 1 OMVARI ANCE . V « 

•  ) IX.3HRMS. 13X.5HRMS/U./.4X.7FI6.7) 

1210  FORMAT! lH0,69X.SHA0JUSTED<9X.BHADJUSTED</< I4X .9HSICMA.9X .6HL0G  IN  )  . 

•  12X  <9HN*Z/U< 7X«SHV*N4F IN) . 7X • 1 3MW*N*F IN )/U*42 .4X . 

4  8HV*N4F IN ) • 7X . 1 3HV *N*F IN ) /U442 ) 

1220  F0RMATI1H  <3X<F16<7< 7X.9H-|I^|NITV.SF14.7) 

1230  FORMAT! IH  .3X.7FI6.7I 

1240  FORMAT! IH) .49HTMCRE  IS  NOT  SUFFICIENT  DATA  TO  MAKE  A  POWER  . 

•  I 6M SPECTRA  ANALYSIS) 

NUMSLOCKaO 

10  READ  1 000 .11 CHUSE 111*1*1.71  S  RE AO  1000. I  I ANUSE 1 1  I  . I ■! . 7  1 

read  1 0] fa. i alpha i i i .beta i i  )« ai  i  i.bi  i  i.ci  i  xoi  i  > <ei  i  ><fi  i  > ,gi  i  ) . 

•  HiD.lal.Tl 
DO  30  1 4 1, 7 

IF  I  I CHUSE 111)  30.30.20 

20  RE AO  1 020. | CMT I  I ) . I WMT I  I  1 • I ANO I  I  I  *  I TNO III. I TCNOI I  I »FZ 111 
30  CONTINUE 

CALL  welcorel 

40  READ  1030. ICROPIOEN! I  I . I ■ I • 10 ) 

READ  lOOO.NBS.AWS.NRF .NPF , M.NSTACC .NSTACECO. INOEXO 
READ  I 040. TFAC . DEL TAT 

CALL  ENOP01NT  A  CALL  YSCTUP 

CALL  PHASESUM  A  |F  INOATA)  270.270.90 

90  00  260  J4|,7 

IF  I  I CHUSE I J I  1  260.260.60 
60  |F  ITFAC-O.OOOOOOI |  TO. TO. BO 
70  GO  TO  190  0  00).  INOEXO 
eo  GO  TO  I  1 10. 120), INOEXO 

90  PRINT  1090  A  GO  TO  130 

100  PRINT  1060  A  GO  TO  130 

110  PRINT  1070. TFAC  A  GO  TO  130 

120  PRINT  IOBw.TFAC 
130  IF  INSTAGE-))  140.140,190 

140  PERaOELTAT  A  T AU"NP*OELT AT 

GO  TO  160 

190  I|4|m/2»44(NSTAGE-) )  A  PER* I l*OELT AT 

T AU* I  I 4NP40EL  TAT 

160  PRINT  1090. ICROPIOEN! I ).!■). 10) 

PRINT  1100 

PRINT  1110. 1CHT I J). IWMTl Jl.FZl J). I ANO I J )« | TNO I J ) • I TCNO I J ) ,M , 

•  OCLTAT .TAU.PER 

NJ4J  A  CALL  OSUBKS 

CALL  LSUBM  A  IF  ITFAC-O. 000000) 1  ISO. 180. 170 


PRINT 

PRINT 


NJ4J 

CALL  LSUBM 
170  CALL  VSUBH 
ISO  CALL  USUBH 

SPEED-VELI J>/NP 
CSS-SS- VEL I J ) 442/NP 
RMS*SORTFIVAR> 


CALL  OSUBKS 

IF  ITFAC-O. 000000) 1  180.180.170 


call  tanon 
SS*VELS<  J ) 

VAR* I SS- VEL I J ) **Z/NP ) /NP 
RMSR*RMS/SPEED 


{ 


iWASfilte 


121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 
S  35 

136 

137 

138 

139 

140 
>41 
142 
144 
144 
1  45 

146 

147 

148 

149 

150 

151 

152 

153 
1*4 

155 

156 

157 
1«« 

159 

160 
161 
162 

163 

164 
16* 
1  66 
167 
1  64 
160 

170 

171 

172 
)  73 
1  74 
1  75 
176 
1  77 
178 
1  76 
180 


PRINT  1 200, SUMl. SPEED. SS. CSS. VAR, RMS. RMSR 

NIC •  M  ♦  1  *  IF  (TF4C-0. 0000001 1  t90. 190.210 

190  PRINT  1120 

00  200  1  ■  1  .NIC 

PRINT  1 ! 30, StGt 1 ) .OK ( I  ) .0KK( 1 1 , LH( 1 ) • UH ( I  I 
200  CONTINUE 
GO  TO  230 
210  PRINT  1 140 

DO  220  1 ■ 1 ,NK 

PRINT  1 150.SIG1 I  1 .OKI  I  1 . QKK ( I  1 ,LH< 1 » . VM C I  ).UH< | ) 

220  CONTINUE 
230  PRINT  1 160 

OO  240  t»| ,NK 

PRINT  I I70.S161 I  ) .FREQ  1 1  I ,FNO< 1 l.FNI I  I , NFN ( t I , NFNR < I  I ,NPFN( I  1 

240  CONTINUE 
PRINT  1180 
OO  250  I • 1 ,NK 

PRINT  l I90.SIGI I  ) .FREQ ( 1  ) , FFN (  1 t .NFFN1 I  ) .NFFNRl I  I .NPFFNI  I  ) 

25 C  CONTINUE 

PRINT  1210 

OO  255  | » I , NX 

IF  < 1  —  t  >  252.252.241 

241  AA»L0GF<FRF0< I ) I 

242  BB*FOEO (11*1 KMT ( J I /SPEED 

CC»VAR*NFN< | )  *  DO*VAR*NFFN( | ) 

CF »RMSR**2*NFN ( | )  *  FF»RMSR4*2*NFFN( | ) 

IF  < I -I  )  253.253.254 

253  PRINT  1 22 0 , S I G ( 1 1. 68 . CC , EE  « 00 • FF 
GO  TO  255 

254  PRINT  1230, S1G( 1 I. AA.Be.CC.EE. OO.FF 
244  CONTINUE 

260  CONTINUE 
GO  TO  280 
270  PRINT  1240 

280  RE AO  1000. 1CR0PCMG  *  IF  (ICROPCHG)  300.300.290 

290  GO  TO  <40. 1 0) , ICROPCHG 
300  FNO 

Subroutine  oatafac 

COMMON/NUMBERS/N8S.NPS.NBF.NPF .m.np.nstage. INOEXQ.NJ.NDATA.TFAC. 

*  oeltat 

COMMON/CHAAWELS/ 1CHUSE<7 ) . I ANUSE  < 7  I . I  SAVE  <7 ) 

COMMON/XSUMS/XX  <31.7|.X1<31.7),XK<3I.7).XS<I 000 . 7 ) . DAT A  < 62 . 7  )  . 

*  XT  <31 .7) .XJ<  3|  .7) ,T J<31  .7 ) .XTBAR<7I .S  <71 
00  20  J«1 , 7 

IF  <  | CHUSE  <  J  >  >  2C. 20. |0 
|C  XTBAP<J)»<XT< Jl-XJI J]*TJ< J)/NP)/NP 
20  CONTINUE 

FNO 

SUBROUTINE  OATASFT< 1 , IPS. JP.T ) 

COMMON/NUMBERS /NBS.NPS.NSF .NPF .M.NP.NSTAGE . 1 NOEXO .NJ.NDATA . TFAC  . 

*  OELTAT 

COMMON/CHANNELS/ | CHUSE  <7 ) , 1 ANUSE  <7 ) . 1SAVE  <  7) 

COMMON/XSUMS/XX  <31. 71. X|  <31. 7), XX <31  .7 ) .XS  < 1 000,7 ) . OAT A  < 62 • 7  )  , 

*  XT<3I.71.XJ<31,7|.TJ<3|.7).XTBAR<7).S<7) 

COMMON/ ANEMF  ACS/U  <  255 • 7 ) .ALPHA  <7 1 »BETA<7).A<7),B<7).C<7).0<7). 

*  E<  7 1 .F  <7 ) .G<  7 1 ,M<7> 

COMMON/FXTRAS/VV  <  7 ) . VEL  < 7 ) . VELS <  7 ) .NSTAGECO 

COMMON/ 1 /N'JMPER.NUMBLOCK .NUMPNTS. I  flag*  1 CHANNEL  < 1680) . I  OAT A < 1680  ) 

I F  <nSTAGF-1 )  2.2.1 
1  NPHASF»<M/2)**<NSTAGF-1 ) 
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181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 
19? 

193 

194 

195 

196 

197 

198 
1  99 
200 
20] 
202 

203 

204 

205 

206 
207 
206 

209 

210 
21  1 
21? 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
218 

239 

240 


10 

1  1 
12 


no  50  j»i .7 
IPS* IPS* 1 
X* I  284 1 0474 ( IPS) 
IF  (NSTAGF-1 1 


*  IF  (ICHUSE(J)I 

*  L* 1 284 1  SAVE  (  J 1 


50.50.10 


13.13.1 1 

VV (  J  1  ■  V  V  (  J  1  +U  ( K  .  J  1 
VEL I J 1  * VEL ( J 1  ♦ VV I J 1 /NPM A  SE 
VELS ( J 1 *VELS ( J 1 4  < VV ( J 1/NPMASE 1**2 
VV(J]*0.0  % 

13  VEL  <  J I  * VEL ( J 1 4U ( X  »  J 1  * 

14  DATA ( I • J1*U(X. J1-TFAC*U(L.J]  * 

20  XT! J]*XT( JI+OATAf 1 . JI*T  * 

TJ(  j]*7j:  J14T  % 

70  S( J]*S( Jl4DATA( 1 ,J1 
40  1  SAVE ( J  1  ■ I OAT  A (  IPS ) 

50  CONTINUE 
FNO 

SUBPOUT INE  OATASUM1 I LOOP .KLOOP .LCMECK I 

COMMON/NUMBER  S/NBf' .  NPS  . NBF . NPF . H.NP.NST AGE • I  NOE  XO .NJ.NOATA.TFAC. 
*  OEITAT 

CONMON/CMANNELS/ 1 CMUSE ( 7 1 . I  AMUSE (71 . I  SAVE (71 

C0MMON/XSUMS/XX(31 .71 .X| (31 .7) .XX (31  .71 .XS( 1 000.71 . OAT A( 62.71  • 


*  IF  (JP-NPKASE)  14.12.12 


GO  TO  14 

VELS ( J  >  *VELS ( J 1  ♦ U (X  «  J  > **2 
IF  (NSTAGE-t 1  20.20.30 
XJ( J1*XJ( J>40*TA( I • J 1 
GO  TO  40 


10 

20 


30 

40 


30 

60 

70 

80 
90 
1  00 


no 

IF 

IF 

00 


XT (31 .71 *XJ(31 .71 .TJ(31 .7] .XT8AP(71 .S(71 

I  00  J*| .7 

( I CMUSE (J) 1  100.100.10 


*  A  J*0 . 5* ( I  I  “  I  1 

*  XT ( J 1 *XT ( J 1 4XS ( N. J 1 *T 

*  TJ( J)*T  .  JI4T 


*  IF  (NSTAGE-I  1  50.50.60 

*  XS  (L  .  J 1  *OAT  A  (L  •  J  1 


»  XX (1.J)*XK(t.J)+XS(L«J1 


(NSTAGE-1 1  40. 40. 20 
30  N* 1 .XLOOP 
1  I  •  ( M/2  1 ** (NST  AGE” 1  1 
T*( I l*N-AJ1*0ELTAT 
XJ(J)*XJ( J1+XS(N.J> 

CONT INUE 
DO  90  1*1. I LOOP 

no  no  x * i .kloop 

L*X«1-1 

XS (X • J 1 *DAT  A (X , J | 

IF  (L-LCMECK)  70.70.90 
XX ( I • J]*XXi 1 . J I 4XS ( X . J I *XS (L • J ) 

XI  (  1  .J)*XI ( 1 . J|4XS(K.JI 
CONT INUE 
CONTINUE 
CONTINUE 
FNO 

SUBROUTINE  ENDPOINT 

COMMON/NUMBEPS/NBS.NPS.NBF.NPF .M.NP.NSTAGE. INOEXO.NJ.NOATA.TFAC. 
*  deltat 

COMMON/EXTR AS/VV (71. VEL ( 7 1 .VELS(7I .NSTAGECO 
NP*(NBF-NBS4]  1  *  1  680-  ( NPS”  I  )-(  1  680-NPF  1 
|F  (NSTAGFCO-I )  1.1.2 

1  NSET»NP/(7*M )  •  NT£MP*NSET*7*M 

GO  TO  3 

2  I  I ■ ( M/2  1 ** (NST AGECO— 1  1  *  NSET-NP/ ( 7* I  I  1 

NTEMP*NSET*7*| | 

3  NREMOVE*NP-NTEMP  *  IF  (NREMOVE— NPF I  20.10.10 

1 0  NBF*NBF-I 

20  NP«NTEMP 

MPF*NP-(N8F-NBS4 I  I  * 1 6604 (NPS- 1  141680 

NP*NP/7 

FNO 

SU5ROUT | NF  FA NON 

TYPE  REAL  lh.nfn.nffn.nfnr.nffnr.npfn.npffn 

COM50N/NUM8ERS/NBS. NPS. NBF.NPF .M.NP.NSTAGE. 1N0CX0.NJ. NDAT A . TFAC . 
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DT ■ I  I  *DELT  AT 


FREQ( I )»K/(2.0*M*0T) 


•  GO  TO  (80.90) . INOExQ 
A  GO  TO  )00 

•  B*P  l*FREQ( I )*TAu 
»  D* (SINF (B)/B)**2 

•  NFN( 1 )«FREQ< | )*FN< I t 

*  A*CUBERTF (FREQ! I ) )**5 

*  NPFFN(  I  )*A*FFn(  I  ) 

*  FFSAvE*NFFn( I ) 


160. 160. ISO 


*  DELTAT 

COMMON/VALUES/OK  <  31  )  .0XX(31  )  .  L  H  (  3  1  ).  UH  (31)  ,FNQ(  31  )  .FN(31  )  .FFNOl  )  . 

*  NFN  ( 3 1  ) . NFFN  ( 3 1  ) .NFNR ( 3  I  ).NFFNR(3I  ).NPfN(3l  )• 

*  NPFFN ( 3 1 ) .VH(31 I.FREOOI ).SIG(3I ) . SUML 

NK«M+1  A  Pt«3. 1415927 

IF  (NSTAGF-l )  10, 10,20 

10  DT»OELT  AT  A  TAU*NP*DELT AT 

GO  TO  30 

20  I | ■ (M/2 )** (NSTAGE-1 )  A  DT*I1*DELTAT 

T  AU* 1 l»NP*DELTAT 
30  DO  160  I»1 .NX 

*•1-1  *  FREO( I  )»K/(2.0*M*0T) 

S|G( I  )»X*OT 
60  DN»1 ,0/(2.0*M#OT ) 

70  FNQ( I )*UH( I )/0N  A  GO  TO  ( 80 . 90 ) . I NOEXQ 

80  FN( 1 )*FNO( | )/QK( 1 )  A  GO  TO  ) 00 

9C  FN( I l»FNO( | l/OKX ( ) ) 

100  AbP|4FR£Q( | )*DT  A  B«P I *FREQ ( I )*T AU 

C* (SINF ( 4 )/A )**2  A  D» (SINF (B)/B)**2 

FFN( | )»FN( 1 )/(C-D>  *  NFN( 1  )»FREQ( I  >*FN( I  ) 

NFFN(  I  )*FR£Q< I  )*FFN(  1  )  A  A*CUBERTF ( FREQ t I ) >**5 

NPFN( I I  ■ A«FN ( I )  A  NPFFN( I )«A*FFN( I ) 

IF  ( 1-1 )  110,1)0,120 

) 1 0  FSAVE-NFN(I)  A  FFSAvE»NFFN( I ) 

GO  TO  160 

120  IF  (NFN( I 1-FSAVE )  140.140.130 
1 30  FSAvE-NFN( I ) 

140  IF  (NFFN( I l-FFSAVE )  160,160.150 
150  FFSAVE«NFFN( I ) 

) 60  CONT INUE 

DO  170  1*1. NX 

NFNR( 1 )*NFN( | l/FSAVE  A  NFFnR( I >»NFFN( I l/FFSAVE 

170  CONTINUE 
END 

SUBROUTINE  INITIAL 

COMMON/NUMBERS/NBS.NPS.NBF .NPF .M.NP.NSTAGE . INOEXQ.NJ.NDATA.TFAC. 

*  DELTAT 

CO*»*»ON/CHANNELS/ICMUSE(7).  I  ANUSE  (  7  )  .  ISAvE(7) 

C0MM0N/xSUMi/XX( 3I.7).X|(3).7).XK(3).7),XS( ) 000 . 7 ) . DAT A ( 62 • 7  )  . 

*  XT (31,7),XJ(31 ,7)«TJ(3I .7) ,XTBAR(7 ) . S( 7) 
COMMON/EXTRAS/VV (7) «  VEL ( 7 ) . VELS ( 7 ) .NSTA&ECO 

N'<«*4  ) 
no  50  j«i .7 
F  (ICHUSriJI)  50.50.10 

10  XT ( J ) *0 • 0  A  XJ(J)*0.0 

T  J ( J I  *  0 • 0  A  VEL (J 1*0.0 

VEL^I J'«O,0 
IF  (NSTAGE-1 )  30.30.20 

?0  s< JI.0.0  A  VV(J)*0.0 

30  DO  40  |* ) .NX 

XX(|.J)«0.0  A  X|(I.J)*0.0 

XX ( | . J)«C.O 
40  CONTINUE 
50  CONTINUE 
END 

SUBROUTINE  LSU"H 

TYPE  REAL  LH, NFN. NFFN. NFNR. NFFNR.NPFN. NPFFN 

COMPON/NUNBERS/'NBS.NPS.NBF .NPF .M .NP.NST AGE . I NDEXQ . N J . NO*  T  A . TFAC . 

*  DFLT  AT 

CCMMON/VALUES/QX (3) I ,QXX( 31 ).LH( 31 ) ,UH(3I ) ,FNQ( 3) ) .FN( 31 ) ,FFN( 3) ) , 


A  NFFNR( I )*NFFN( I l/FFSAVE 
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30 1 

302 

303 

304 

305 

3Cft 
307 
300 
30<? 
3)0 
31  1 

312 

313 

314 
31*< 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 
34  0 
34  1 
342 
3*3 

344 

345 

346 
34  7 

348 

349 

350 

351 

352 
333 

354 

355 
336 
357 
338 
739 
760 


'  NFNI31  )  ,NFFN(3 I  )  .NFNROl  )  ,NFFNfl(3l  ).NPFN(3I  )  . 

•  NPFFN (31 ) « VH (31  ).FRE0<3!  )  «SIG(3I  ) . SUML 

DIMENSION  0(31) 

NK«M* 1  «  NM-NK-1 

P I ■ 3, • 4  1 5927  4  SUML a  0 , 0 

DO  JO  1*1. NX 
GO  TO  ( 10«20 ) ( INDEXO 

)0  0(1 ) *OK (  |  )  *  60  TO  30 

20  0( I )  *QKX  (  |  ) 

30  CONTINUE 

no  no  i«i(Nr 

IF  ( 1-1 )  50.50(40 
40  IF  (l-NKI  70.90,90 

90  LH( 1 )*(0( 1 )*0 (NX ) )/(2.0*M)  S  SUMaQ.O 

no  60  J*2 ,NM 
SUM* SUM4Q ( J ) 

60  CONTINUE 

LM( | )*LH( I )*SUM/M  S  SUML«SUML*LH( I ) 


GO  TO  1 1 0 

70  LH ( I ) a ( 0 ( 1  no (NX )*C0SF (  ( 1-1  )  *P  I  )  )/M 
‘UMaO.O 


DO  80  J  a  2 . NM 

AN6a(  |-|  )•(  J-t  )40|/M  *  SUM*5UM«.Q(  J1*C0SF(ANG) 

80  CONTINUE 

LH(  I  )a|_H(  |  )>2,0*SUM/M  •  SUMLaSUML4LH(  I  > 

GO  TO  11 0 

90  LM(  1  )«(0(  )  (♦(-)  )**M*0(NX )  )/(2.0<>M) 

SUMaO.O 

DO  100  Jap.NM 

SUMa  SUM  4-  ( — )  >**(J-1  1*0  (  J  ) 

100  CONTINUE 

LH(  I  )*LH(  1  nSUM/M  •  SUMLaSUML*LH(  I  1 

I  10  continue 

END 

SUBROUTINE  PHASE Sum 

C0MM0N/NUMBERS/NBS,NPS.N8F,NPF .M.NP.NSTAGE, INOExO.NJ.NDATA.TFAC. 

*  DELTAT 

COMMON/ CHANNEL S/ I CHUSE ( 7  > , I ANUSE ( 7 ) , I S A VE ( 7 ) 

COMMON/XSUMS/XX (3)«7),X| (31 ,7) ,XK (31 , 7 ) , XS ( 1 000 . 7 ) . DATA ( 62 , 7 ) , 

•  XT(3> ,7) ,XJ(31 ,7) ,TJ(31 ,7) ,XTBAB(7) ,S(7) 

COMHON/|/NU**PEB.NUM5LOCK  .NUMPNTS,  I  FLAG,  ICHANtCL  (  1680)  •  I  DATA  (  1680  ) 
NX |  •  IF  (NSTAGE-1  >  10,  If/, 20 

)0  NPHASEaM  •  GO  TO  30 


?0  NPHASEa ( M/2 ) ** ( NSTAGE - 1 ) 
30  CALL  INITIAL 
*0  NT 1 ME»N8S— NUM8L0 CX 
DO  50  | a) , NT 1  ME 
CALL  CRACX 
30  CONTINUE 
60  IPSaNPS-1 
70  IPF* | 680 
80  IPFaNPF 
90  NSa| 

JP»0 
I  00  NPaO 


*  IF  (NUM8L0CX— NBS )  40,60.60 


•  IF  ( NUM8L0CX-NF/F )  70.80.80 

•  GO  TO  90 

•  NFa2*NK 

•  IF  (NSTAGE-I >  110.110.100 


1 |0  00  260  laNS.NF 


jpajp«| 

120  NSTORE* I 
60  TO  140 
170  TaO.O 


•  IF  (NSTAGE-1)  120,120.130 

•  TaJPODELTAT 
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i 

i 


m 


361 

140 

362 

150 

363 

1*0 

364 

368 

366 

1  70 

3*7 

1  80 

368 

1  90 

369 

200 

370 

210 

371 

220 

372 

230 

373 

374 

240 

375 

280 

376 

260 

777 

778 

270 

779 

280 

380 

781 

290 

782 

383 

300 

384 

385 

310 

386 

320 

387 

388 

389 

330 

790 

39| 

340 

792 

780 

393 

394 

3*0 

395 

79* 

370 

397 

398 

380 

399 

390 

400 

400 

40  | 

410 

402 

420 

403 

404 

4O5 

40* 

4O7 

408 

409 

410 

41  1 

412 

41  3 

4  1  4 

415 

416 

10 

4  )  7 

418 

20 

4  1  9 

30 

42^ 

40 

CALL  DATASET ( | , IPS. JP. T ) 

IF  (JP-NPHASEI  190. I *0.160 
NP*NP*I 
00  180  L ■ I  .7 

IF  (ICMUSE(L))  180.180,170 
XSINP.L  >»S(L l/MPHASE 
CONT INUE 

IF  (IPS-IPFI  2*0,200.200 
IF  (NUMBL0CF-N8F )  230.210,210 
IF  (NSTAGE-I I  220.220,290 
nST0P»2 
CALL  COACtC 

IF  (NUM8L0CF-NBF |  240.250.250 
I PF  » | *8  0 
|PF«NPF 
CONTINUE 
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